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ABSTRACT 

We study the evolution of the orientation of the orbital plane of massive black hole 
binaries (BHBs) in rotating stellar systems in which the total angular momentum of 
the stellar cusp is misaligned with respect to that of the binary. We compare results 
from direct summation A^-body simulations with predictions from a simple theoreti- 
cal model. We find that the same encounters between cusp stars and the BHB that 
are responsible for the hardening and eccentricity evolution of the binary, lead to a 
reorientation of the binary orbital plane. In particular, binaries whose angular mo- 
mentum is initially misaligned with respect to that of the stellar cusp tend to realign 
their orbital planes with the angular momentum of the cusp on a timescale of a few 
hardening times. This is due to angular momentum exchange between stars and the 
BHB during close encounters, and may have important implications for the relative 
orientation of host galaxies and radio jets. 
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1 INTRODUCTION 

In the standard ACDM cosmological scenario, massive black 
hole binaries (BHB s) are the natural out come of mergers 
of massive galaxies fSegelman et al.ill980l ). The orbital de- 
cay of BHBs has been studied in detail under the assump- 
tion of purely stella r systems (e.g. Makino fc Funatol 20041: 
Berczik ct al. 2005; Baumgardt et ajJ l2006l : iMerritt et al.l 
200 7 ; Scsana ct al. 2007; Scsana .20101) as w e ll as g aseous 



environment s (e.g . lArmitage fc Nataraian 1200! 



Escala et al.l 



2004, 12005 



Cuadra et al.l 



Dotti et al.l l2006l. 1200' 



2005 



200S 



2003 : iLodato et allbOOalRoedig etallboill ) 



Most studies of the evolution of BHBs in stellar envi- 
ronments consider non rotating, spherically symmetric nu- 
clei, despite t he natural occurrence of rot ation in merger 
remnants (e.g. iGualandris fc Merrittll201ll ). Simulations of 
MBHs in spherical stellar systems predict a stalling of the 
binary evoluti on at parsec scale separations , the "final par- 
sec problem" (|Milosavlievic fc Merrittlliooal ). which may or 
may not be avoided in the presence of gas dissipation. A 
BHB can, however, reach coalescence even in a completely 
gas poor environment, if embedded in a non-spherically 
symmetric stellar distribution. Purely dissipationless simu- 



lations of non-spheric al galaxy models (iBerczik et al 



iBerentzen et al.ll2009l ) or galaxy mergers ( Khan et al 
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iPreto et aLll201ll : [Cualandris fc Merrittll201ll ) find continu- 
ing hardening of the binary down to separations where en- 
ergy loss due t o emission of gravity waves be comes efficient. 
Simulations bv IGualandris fc MerrittI (|201ll ) show that effi- 
cient hardening is achieved even in systems close to axisym- 
metry. This departure from spherical symmetry is due to 
the rotation introduced by the merger, which results in an 
enhanced rate of stars interacting with the binary. 

A comprehensive and systematic study of the evo- 
lution of BHBs in rotating stellar cusps is still missing. 
ISesana et al.l l|201ll) (hereafter paper I) studied the evolu- 
tion of the eccentricity of unequal mass BHBs in rotating 
systems by means of an hybrid analytical/3-body scattering 
formalism as well as full A''-body simulations. They found a 
strong dependence of the eccentricity evolution of the binary 
on the degree of co-rotation in the stellar cusp, with binaries 
increasing their eccentricity in cusps containing a significant 
fraction of counter-rotating stars. This is relevant to the case 
of galactic minor mergers, in which counter-rotating systems 
can be produced as a result of rotation in the larger galaxy, 
thereby forming very eccentric binaries. 

Another important and potentially observable param- 
eter is the orientation of the binary orbital plane. Its evo- 
lution has been studied in spherical and isotropic systems 
(|Merrittll2003: IGualandris fc MerrittI [20071) . Due to the lack 
of any preferential direction in these systems, the orbital 
plane of the binary can only undergo a random walk, result- 
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ing in small changes in the orbital plane on long timescales. 
Here we study the evolution of the orbital plane of BHBs 
in rotating stellar systems, where, as will be discussed, a 
significantly larger reorientation is expected. 

The letter is organized as follows. In section[2] we de- 
scribe the evolution of the binary orbital plane in spherically 
symmetric models. In section[3] we describe the case of ro- 
tating models, comparing results from A''-body simulations 
with a simple analytical model. Final remarks are presented 
in Section|4] 




2 NON ROTATING CUSPS 

The orientation of the orbital plane of BHBs embedded in 
non rotating, isotropic stellar nuclei has been studied for 
the first time by Merritti (^2002)- In such systems, repeated 
encounters with passing stars cause the orientation of the 
binary to u ndergo a random walk (a "rotational Brownian 
motion", in iMerritd (|2002t )l. From simple analytical argu- 
ments, the expected change in orientation on one hardening 
time Th = |a/a| of the binary (wh ich is of the order of 
120 i nitial binary orbital periods) is (j Gualandris fc Merritti 
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where q = M2/M1 is the binary mass ratio, Mbhb the total 
mass, and e the orbital eccentricity. Numerical scattering ex- 
periments have constrained the normalization in Eq. [1] only 
for circular binaries. Neither the dependence on the eccen- 
tricity nor on the mass ratio has been tested numerically. In 
the case of a binary with AIbhb = lO^m, and q = 10~^, the 
expected reorientation is A6 ~ 9°. 

Here we check the validity of Eq. \T\ by performing A^- 
body simulations of binaries with different eccentricities 
(in the range — 0.99) embedded in non-rotating stellar 
cusps. We use the d irect summation A^-body code <^GRAPE 
jHarfst et al.l |2007|) in com bination with the Sapporo li- 
brarv (jCaburov et ahllioogl ) to accelerate the computations 
on GPU hardware. In our integrations we adopted Mi = 
10'' M0 and q = 1/81 and set the initial semi-major axis to 
at = 0.06 pc. The MBHs are embedded into a stellar cusp 



following a Bahcall-Wolf p(r) 



-7/4 



density profile at dis- 



tances smaller than 1 pc, with total mass Mc ~ 2.5 X lO'' Mo 
and a mass enclosed in the binary orbit equal to 2M2. The 
cusp is modelled with A'' = 32k equal-mass particles, result- 
ing in a black hole to star mass ratio of m*/A/i — 7.5 x 10~^. 
Fig [1] shows the degree of re-alignment of the binary plane 
as a function of the initial eccentricity. The dashed line in- 
dicates the dependence expected from the analytical model, 
while the points show the results of our A'^-body integrations. 
The agreement is remarkable. 

For a given binary, the re-orientation in a non-rotating 
system does not depend only on the mass density and ve- 
locity field of the stars, but also on the mass ratio between 
the binary and a single star in the cuspQ. A binary embed- 



^ Note that this ratio is properly defined only for a single mass 
stellar cusp. In the presence of a mass spectrum, the mean mass 
of the interacting stars is a good proxy for m* . 



Figure 1. Change over one hardening time in the orientation 
of the angular momentum of the binary, as a function of the 
initial eccentricity of the BHB. The dashed line represents the 
theoretical dependence given in Eq. [T] normalized to match the 
data at e = 0.25. The amplitude of the reorientati on in our runs i s 
however comparable to the analytical estimates in lMerrittlll2002l) . 



ded in a stellar cusp would, on average, experience a greater 
realignment as a result of few encounters with massive stars 
than as a result of many encounters with low-mass stars. 
This feature is typical of any random walk process, in which 
there is no preferential direction for the change of the binary 
plane in a single encounter and changes due to different en- 
counters tend to cancel each other out. The same random 
walk behaviour is observable in rotating cusps, if the angu- 
lar momentum of the binar y is aligned with that of the cusp 
(|Amaro-Seoane et al.ll2O10l ). This behaviour is illustrated in 
the lower panel of Fig |2l which shows the dependence of 
the change in 6 on the number of particles used to model 
the cusp, i.e., on the rrit / Mbhb mass ratio. The agreement 
between the numerical results (full circles) and Eq. [T] (solid 
line) is outstanding. 



3 ROTATING CUSPS 

Here we study in detail the re-orientation of a massive BHB 
embedded in a cusp with net rotation. To this aim, we per- 
form high accuracy A'-body simulations similar to those pre- 
sented in the previous section, the only difference being the 
degree of rotation that we enforce in the cusp. This is ob- 
tained, as in Paper I, by reversing the sign of all velocity 
components for a random subset of cusp stars, at the time 
where the BHB is added. In particular, we generated mod- 
els with a fraction = 0.875 of co-rotating stars and varied 
the initial angle 9o between the total angular momentum 
of the cusp and the angular momentum of the binary. We 
considered cases with 60 = 0°, 30°, 60°, 90°, 120°, 150°. 

The evolution of 6 as a function of time is shown in 
Fig. [3] (solid lines) for all runs. Binaries whose angular mo- 
mentum is initially mis-aligned with respect to that of the 
cusp tend to realign it on a time-scale of a few hardening 
times. The re-alignment can be quite significant, with an 
average change in 9 of the order of 50% of the initial value, 
and up to 70% in some cases. For example, the case with 
9o = 150° reaches a final value of 6^ ~ 45°, with a change 
of more than 100° over one hardening time. This is about 
100 times larger than the reorientation measured for a sim- 
ilar isotropic model (see Fig. [1]). As will be discussed be- 
low, such large reorientation of the orbital plane may have 
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Figure 2. Change in the binary plane orientation per hardening 
time versus number of particles in the cusp (in the range N = 2h 
to N = 128k) or, equivalently, ratio between star mass and MBH 
mass. In all the runs e = 0.5. The lower panel refers to isotropic 
models, while the upper panel refers to rotating models with = 
0.875 and an initial angle between the binary angular momentum 
and that of the stellar cusp = 30°, 90°, 150°. 



Figure 4. Solid lines: Distribution of the angle dcusp between 
the angular momentum vector of a single star and the initial 
angular momentum vector of the cusp, for the run with Oq = 150° . 
Only unbound stars are considered. The different curves refer to 
different times during the evolution, from 1 Ti; to 8 T^ . Stars are 
ejected roughly isotropically from early times. Dotted line: Initial 
distribution of (9cusp for all cusp stars, normalized by a factor of 
10, for an easy comparison. 
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Figure 3. Evolution of the angle 9 between the angular momen- 
tum vector of the binary and that of the stellar cusp. Time is 
expressed in units of the initial binary orbital period (bottom) 
and of the hardening time (top) . Different solid lines are for mod- 
els with different initial values of 0, from Bq = 0° to do = 150°. 
Dotted lines represent the predictions from the theoretical model 
while dashed lines indicate the predictions modified to take into 
account the finite reservoir of angular momentum in the stellar 
cusp. 



observational consequences. However, we note that the re- 
alignment is not complete but stops at a finite value of 9. 

The evolution of 6 can be explained in terms of a simple 
analytical model. Let Ib = Lb / l-ihe the angular momentum 
per unit mass of the binary, with ^ — A4iM2/{Mi +M2) the 
reduced mass, and — Lt/m* the angular momentum per 
unit mass of a star interacting with the binary. Conservation 
of angular momentum during an encounter (assuming m* <C 
Mi+ M2) gives 



AIb = —AL 



(2) 



Considering that, once normalized per unit mass, the 
average change per encounter in the stellar angular momen- 
tum i s of the order of the binary angular momentum (|Merrittl 
|2002| 1^. each interaction contributes, on average, a factor 



< AIb >~ — Ib 



(3) 



to the binary angular momentum. The rate at which Ib 
varies with time is 



dlB 

dt 



dric 



dt 



< AIb > , 



(4) 



where dricnc/dt is the rate of binary/star interactions. The 



^ This allows us to ignore the evolution of the magnitude of the 
binary angular momentum due to hardening when modelling the 
evolution of the plane orientation. The change in 6 would be the 
same even if wo considered the decay of the binary orbit. 
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Figure 5. Sketch of the model explaining the re-alignment of the 
binary plane. Lb o the initial angular momentum of the binary 
(perpendicular to the orbital plane, shown for clarity in the left 
panel). i is the angular momentum of the binary after re- 
alignment. Ai is the amount of angular momentum transferred 
at any interaction with a single star (dashed line). Note that we 
assume parallel to the total angular momentum Lcwsp 

of 

the cusp. 9q and Si indicate, respectively, the initial and final 
angles between i/g and J^cusp. 

Panel (a) shows a case with iB,o 
partially aligned with icuap (^0 < 90°). Panels (b) and (c) refer 
to the opposite case {Oq > 90°). Vectors are not on a physical 
scale. 



rate can be estimated as 



d72cnc j-j2 

= TT Tliiti 2 Wrcl 
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where n, is the number density of stars in the cusp, i?i,2 is 
the influence radius of the secondary hole, and Urei is the 
mean relative velocity between the stars and M2. Note that 
we have implicitly assumed that: i) only stars efficiently in- 
teracting with the secondary modify the angular momentum 
of the binary; ii) the dynamics outside the influence radius 
of the secondary is completely dominated by the primary. 

In order to apply Eas. l2l5l we need to specify the aver- 
age direction of < AZb >• This is opposite to the direction 
of the variation of the angular momentum of the star during 
the interaction. In a reference frame in which the z direction 
is aligned with the total angular momentum of the cusp, the 
average x and y components of the angular momentum of 
the stars before an encounter are zero. In this simple model 
we assume that the final angular momentum of a star that 
experienced an interaction is isotropically distributed. This 
is justified by the observation that, at any time, stars after 
an interaction with the binary are much more isotropically 
distributed than before, as shown in Fig. |4] for the popu- 
lation of unbound stars. As a consequence, the variation 
< Al, > — < Z,_i — l^fi > points in the direction opposite 
to the angular momentum of the cusp Zcusp, and < AIb >, 
averaged over repeated encounters, points in the z direction, 
parallel to icusp- The fast alignment that we observe in the 
A'^-body simulations is a direct consequence of such prefer- 



ential direction of AIb ■ Fig. [5]shows a schematic view of our 
model, for different initial orientations of the binary. 

The evolution of predicted by our model is shown with 
dotted lines in Fig. [3] Naturally, the case with Ib aligned 
with icusp does not show any evolution. Despite the simpli- 
fying assumptions made, the agreement between the initial 
evolution of 9 in the runs and in the model is very good. 

The predictions of the model and the results of the 
Ai'-body runs disagree at late times. In the simulations 
the binary stops re-orienting well before perfect alignment 
is reached, while the model predicts 0^0 asymptot- 
ically. This difference can be explained simply in terms 
of a depletion of stars in t he centre of the cusp. As in 



many previous studies (e.g. Milosavlievic fc Merritt 
Milosavlievic et al.ll2qd3: lMerrittll2006l : iMerritt et al. 



2001 



2007 



Gualandris fc Me^ritt^l2011^ ■ we observe the formation of a 



central core due to slingshot ejections of stars. After a core 
has been excavated in the centre, the rate of encounters dras- 
tically drops and the evolution of the BHB's orbital plane 
halts. This translates into a maximum angular momentum 
^max that can be transferred from the stars to the binary . 
Since in our model we assumed the transferred angular mo- 
mentum to depend on the properties of the cusp and not 
on those of the binary (the stars are ejected isotropically), 
^max must be the same in every run. We computed the value 
of /max that needs to be transferred to the binary to match 
exactly the evolution of the ^0 = 120° run, and we checked if 
such a limiting angular momentum can account for the other 
runs' behaviour. The results of this test are shown in Fig.|3l 
with dashed lines. The agreement with the A'^-body results is 
greatly improved by the introduction of a maximum angular 
momentum that can be transferred to the binary. 

As a final note, we emphasize that the re-alignment in 
rotating cusps does not depend on the m, / Mbhb ratio, as 
opposed to the non-rotating case. Since the rate of encoun- 
ters under every assumption depends on the number density 
of stars, and the amount of angular momentum transferred 
per encounter is proportional to m, , d^B /At depends on the 
mass density p of the cusp, but does not show any other 
dependence on m«/MBHB. This is a consequence of the fact 
that the alignment in each encounter has a preferential direc- 
tion (i.e. towards the angular momentum of the cusp). The 
upper panel in Fig. [5] shows the alignment AO experienced 
over an hardening time by binaries in rotating cusps, as a 
function of the number of stars in the cusp. Different lines 



are for different initial angles between Ib and lc 



While 



in the non-rotating case the ratio between the re-alignments 
in the lowest and highest mass resolution runs is ~ 30, in 
this case the re-alignment shows only statistical fluctuations, 
with no systematic trends. 



4 CONCLUSIONS 

We studied the evolution of the orientation of the orbital 
plane of unequal mass BHBs in isotropic and rotating stel- 
lar systems. The evolution in isotropic systems is character- 
ized by a rand om walk on long time-scales, as predicted by 
iMerrittI (|2002l ). The evolution in systems in which a frac- 
tion of stars is counter- rotating with respe ct to the binary 
was briefly discussed in ISesana et al.l l|201ll ). If the fraction 
of counter-rotating stars is small, the change in the orbital 
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plane is consistent with the change expected for isotropic 
cusps. If, on the other hand, a large fraction of stars is 
counter-rotating, the plane of the binary evolves consider- 
ably and the binary angular momentum even reverses in the 
case of a fully counter-rotating cusp. 

Here, we investigated for the first time the evolution of 
the binary plane in rotating systems in which the binary 
angular momentum is initially misaligned with respect to 
that of the stellar cusp. We find that mis-aligned binaries in 
rotating cusps tend to align their angular momentum with 
that of the cusp on a time-scale of a few hardening times. 
The alignment process is linear with time, hence, unlike the 
random walk observed in fully isotropic systems, does not 
depend on the mass ratio between the MBHs and the stars. 

The realignment of the binary plane that we observe in 
the simulations may have significant implications. The direc- 
tion of the spin axis of the single MBH that results from the 
merger of the holes due to emission of gravitational waves 
is affected by the orientation of the binary plane before co- 
alescence. The spin axis, in turn, determines the orientation 
of the accretion disk aro und the remnant b lack hole via the 
Bardeen-Peterson effect (jBardeen fc PettCTS on 1975) and, in 
radio-loud AGNs, the direction of the radio jet. 

After the coalescence of a BHB, the spin of the rem- 
nant is not necessarily aligned with the spins of the two 
progenitors, and jj. If one of the two MBHs in the bi- 
nary is radio-loud (and the remnant stays radio-loud), the 
change in spin dire ction would result in a reorientation of 
the relativistic je10 (|Merritt fc Ekersll2002l '). The orientation 
of the orbital plane of the binary has implications for the di- 
rection of J J,. The interaction with a rotating cusp tends to 
force the binary to co-rotate with the central cusp. At coa- 
lescence, the binary orbital angular momentum contribution 
to jj, will result in a spin (and possibly a jet) preferentially 
aligned with the cusp angular momentum. 

We do not expect full alignment, however, because: (i) 
the alignment we find in the simulations is not complete, and 
could halt in the absence of further refilling of stars or in the 
presence of isotropic refilling, leaving misalignment angles of 
up to 60°, (ii) in a very unequal-mass binary (e.g. q — 1/81, 
like the o ne we study), depe nds strongly on Using the 
results of lRezzoUa et al.l l|2008h . if ji ~ 1, the angle between 
J J and jj. is of only a few degrees, while ji « 0.1 results in 
an angle of > 20°. 

Such a preferential alignment is in fact obs erved. In 
weak radio-loud AGNs, iBrowne fc Battvd l|2010f) found a 
tendency for the axis of the radio emission to align with the 
minor axis of the starlight of the host, an oblate rotationally 
supported elliptical. By contrast, they found no preferred 
radi o-optical alignment among the ra dio-louder objects, (see 
also lSaripalli fc Subrahmanvanl20p"^'). possibly hosted in tri- 



axial non-rotating ellipticals (e.g. Kormendv et al.ll2009f ) 



^ This is valid in gas poor mergers. In gas rich environments, gas 
accretion causes the spins of the MBHs to align with the angu- 



lar m omentu m of the bina ry (Bogdanovic et al.ll2007l :[ Dotti et al.l 
12010 : Volontcr Tet al.l[2oi O: Kcsdcn ct al.l l2010l '), preventing any 
drastic change in the spins direction at coalescence. 
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